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We investigate low-energy properties of the alternating 
spin chain model composed of spin si and S2 with a singlet 
ground state. After examining the spin-wave spectrum in 
detail, we map low-energy spin excitations to the 0(3) non- 
linear sigma model in order to take into account quantum 
fluctuations. Analyzing the topological term in the resulting 
sigma model, we discuss how the massless or massive exci- 
tations are developed, especially according to the topological 
nature of the alternating spin system. 

PACS: 75.10.Jm, 05.30.-d, 03.65.Sq 



I. INTRODUCTION 

One-dimensional quantum spin systems show various 
interesting phenomena, having been continually provid- 
ing hot topics. The most important example is the spin 
one-half antiferromagnetic Heisenberg chain, which was 
exactly solved by the Bethe-ansatz method. It is now well 
known that this model shows critical behavior specified 
by the level-1 Wess-Zumino-Witten model. Keeping the 
integrability by adding multiple terms, higher-spin mod- 
els were constructed and their critical behavior was 
also classified. Progress in a slightly different direction 
was made by Haldane , who conjectured that general 
higher-spin models behave in quite distinct ways accord- 
ing to whether the spin of the models is integer or half- 
odd integer. The massive phase of integer spin models, 
so-called the Haldane phase, has been indeed observed by 
experiments, and now is understood as the valence-bond- 
solid state The problem whether a given quantum 
spin system is massless or massive provides nontrivial 
and interesting issues, which have been stimulating the 
study of quantum phase transitions. 

More recently, the mixed spin chains with an alter- 
nating array of two kind of spins have attracted much 
current interest 0,^. The systems experimentally found 
so far 10] show a magnetically ordered ground state even 
for the antiferromagnetic case. For such ordered systems, 
quantum fluctuations may not play a crucial role. How- 
ever, it is quite interesting to address the question what 
happens when the magnetic order disappears due to frus- 
trations, etc. If such systems are realized experimentally, 
quantum fluctuations should certainly play an important 



role, exhibiting interesting phenomena such as the com- 
petition between the massless and massive phases. These 
systems may thus provide a new interesting paradigm of 
quantum spin chains. In this connection, an exactly solv- 
able alternating spin chain, which has a singlet ground 
state, has been studied in detail jo). 

The alternating spin models can also be regarded as a 
high concentration limit of magnetic impurities in spin 
chains. In this respect, the problem of the alternating 
spin models is closely related to that for the spin chains 
with impurities, which have been extensively studied for 
e.g. the two-leg ladder system [p^|-p^, the spin-Peierls 
system etc. In fact by changing the concentration 

of periodic array of impurities, we can naturally interpo- 
late the alternating spin systems and spin systems with 
dilute impurities, and discuss their characteristic proper- 
ties in a unified way jist . 

In this paper, we study low-energy properties of the 
mixed spin chain models with an alternating array of 
two kind of spins. We deal with two kind of alternating 
spin models with a singlet ground state. By introduc- 
ing quantum fluctuations to the spin-wave excitations, 
we investigate low-energy properties of the system by re- 
sorting to 0(3) non-linear sigma model techniques. We 
then discuss how the massless or massive excitations are 
developed in the alternating spin models by analyzing 
the topological term in the sigma model. 

The paper is organized as follows. In Section II, two 
types of the models for alternating spin chains are intro- 
duced, and some rigorous statements are given for the 
excited states as well as the ground state. In Section 
HI, we calculate the spin-wave dispersion relation, and 
in Section IV, by converting the model to the non-linear 
sigma model, we study the low-energy characteristics of 
the spin excitations. The last section is devoted to sum- 
mary of the paper. 

A brief report on preliminary results for this paper is 
given in [ p^ , where the effects of massive modes were 
not taken into account in the nonlinear sigma model ap- 
proach. 



II. MODELS AND BASIC PROPERTIES 

In this section, we first introduce two kind of alternat- 
ing spin models with a singlet ground state, and then 
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study them by means of Lieb-Schultz-Mattis (LSM) the- 
orem. 

The models we study in this paper are the quantum 
spin models with nearest neighbor interaction in one di- 
mension. The corresponding Hamiltonian is defined by 

H = ^JjS,-S,+u (2.1) 

i=i 

where Na is the number of lattice sites, and the spin at 
the j-th site is denoted by Sj. We will investigate the 
following two types of alternating spin-chain models by 
setting 



Model A: 



si for j = 1 
S2 for j = 2, 3 
J, = J, 
Na = 67V, 
Model B: 

f si for .7 = 1,2 
•' \ S2 for j = 3, 4 

^(l-7i) 



(mod 3) , 



(2.2) 



(mod 4) , 
for j = 1 



J, = { J[l + (71 + 72)/2] for J = 2, 4 (mod 4), 



J(l-72) 
iV„ = AN. 



for j = 3 



(2.3) 



In both cases, the spin chain is composed by two kind 
of spins si and S2, which can take either integer or half 
integer spin. Note that we also include the effect of the 
bond-alternation for the Model B, which was originally 
introduced by Tonegawa et al. In what follows, we 
assume the antiferromagnetic couplings J > and < 
71-72 < 1- 

Let us start by specifying the ground state properties. 
Since spins are on a bipartite lattice, we can apply the 
Marshall theorem 1 17- 19,1^] to ( |2.l|) , that is, we prove 
that the ground state is singlet without degeneracy for 
both models. Next consider the properties of excitations 
above the ground state by applying the LSM theorem 
|pO| . Note that the Hamiltonian is invariant under 3- (4- 
)site translation for the Model A (B). Therefore, together 
with the fact that the ground state is unique, we can 
construct an excited state of 0{1/N) when Suc = odd- 
integer, where Suc is the summation of the spins in the 
unit cell 



si + 
2(si 



2s2 for A 
-f S2) for B 



(2.4) 



According to the above formula, we can state that 
the Model A can be gapless if si is a half-integer. In 
other cases where si is an integer for the Model A, or si 
and S2 are arbitrary for the Model B, we cannot claim 
anything by the LSM theorem, which suggests that such 
systems may be gapful, as is the case for the ordinary 
Haldane system. 



III. SPIN- WAVE SPECTRUM 

So far we have seen that the Model A can show dif- 
ferent phases, according to whether si is integer or half- 
integer. However, the LSM theorem itself cannot distin- 
guish a gapless phase from a gapful phase with degen- 
erate ground states. In order to clarify this point, we 
wish to use a complementary approach based on effec- 
tive field theory. For this purpose, we first investigate 
the dispersion relation of the spin- wave spectrum for the 
model, which allows us to judge whether the system can 
be mapped to the non-linear sigma model. The next sub- 
section is devoted to deriving the spin- wave Hamiltonians 
by the Holstein-Primakoff mapping, and they are diago- 
nalized in the subsection B. Wc then show some examples 
of the spin- wave spectrum calculated numerically. 



A. Spin- wave Hamiltonian 

1. Model A 

For the usual uniform spin chain (si = S2 case for the 
present model), it is sufficient to introduce the two kind 
of bosons, corresponding to bipartite lattices. However, 
for the model with si ^ S2 we have to introduce the six 
kind of bosons such that 



S, 



= (-)' 



(m) 



s, 



a, 



2,3 

(m)t 
6j+m 



(m) (m) 



2sm - n, 



(m) 
6j+m 



for odd m 
for even m 



(3.1) 



where j = 0, 1, • • • , iV — 1, m = 1, 2, • • • , 6, Sm = si (52) 
for TO = 1,4 (otherwise) and n^"^^ 



6j+m 



(™)t (m) J 



what follows, we approximate the boson Hamiltonian up 
to quadratic order in order to derive the dispersion rela- 
tion of spin wave excitations. In this assumption, the 
resultant Hamiltonian can be diagonalized in the mo- 
mentum space. For this purpose, let us introduce the 
Fourier-transformed operators 



a, 



(m) 
'6j+m 



N~l 



^ exp 



k=0 



27ri k , 
^6 



{6j + m) 



(m) 
'■k ' 



(3.2) 



with N being defined in eq.(2.2). Substituting these 
equations, we find 



N-l 6 



k=0 Lm=l 



+a/A/,„aJ„ + aiAijnam), (3.3) 
where we have used a notation. 
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Am) 



for odd m 
for even m 



(3.4) 



for short, and 6x6 matrices h and A are given by 

1 + a 1 + 



h = Js2diag 1, — , 1, 



f V^p 
\fap p 
J Si p y/ctP 

-Jap yfap 
\fap 

V P 



P 

y 



(3.5) 



We have defined in these equations 



fi 
S2 '' 



(3.6) 



with p = {2ti /Na)k. Matrix elements not written exphc- 
itly in eq.( |3.5| ) are all 0. Note that the operators (3.4) as 
well as the matrix A depend on the momentum, though 
we have not explicitly denoted it. 



The spin-wave dispersion relation for this model is 



uj± = ±\si - S2I + [(si - 52)^ + 4siS2 sin^p] 



1/2 



(3.9) 



Thus the excitation spectrum is given by a quadratic 
function for small momentum. It is only the case for 
si = S2 that it has a linear dispersion. 
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FIG. 1. Spin-wave spectrum of the Model A in unit J as 
functions of the momentum p for si = 1/2 and S2 = 1. 



2. Model B 

In a similar way to the Model A, we now write down the 
spin- wave Hamiltonian for the Model B. Introduce the 
Holstein-Primakoff m app ing with the use of four kinds of 
bosons similar to eq.(3+) by replacing 6-^4. With the 
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same Fourier-transformed operators as eq.(3.2), we have 
the spin- wave Hamiltonian similar to eq.(3^) for which 
h and A are 4x4 matrices defined by 

h = -^^^diag (aFi -I- F, aFi -I- F, aV + F2 , aV + F2) , 



JS2 



/ aVip ^/aTp\ 



aVip \/aTp 
Varp T2P 
\^Tp T2P 



(3.7) 



where F^ = 1 — 7^ and F = 1 + (71 -I- 72)/2. 



B. Diagonalization 

In the previous subsection, we have derived the spin- 
wave Hamiltonians up to the quadratic order in boson 
operators. We can diagonalize these Hamiltonians by 
the Bogoliubov transformation. 

Before discussing the spectrum of the present models, 
let us recall the alternating spin-chains with ferri mag - 
netic ground state |^, i.e, the same model given by (2.1), 
but with 



0.5 




S2j-1 — Si, S2j 



S2, 



J 3 = J- 



(3.8) 



FIG. 2. Spin-wave spectrum of the Model B in unit J as 
functions of the momentum p for Si = 1/2 and S2 = 1. 
Bond-alternation parameters are 71 = 72 = 0. 

Now turn to the present model. We show the spectrum 
obtained numerically in Figs.^ and ||. We can see that 
the dispersion relation ujq oc sinp defined in p G [0, tt] for 
the uniform chain is folded 3 (2) times for the Model A 
(B), and interactions among such modes produce gaps at 
the crossing point. Though we here presented only the 
case with si = 1/2 and S2 = 1, the qualitative features 
are the same for other cases si ^82- Contrary to this, 
in the case si = S2, interactions between the different 
modes disappear and the single dispersion ljq appears. 
What is remarkable is that the lowest mode has a linear 
dispersion, in contrast to the ferrimagnetic case, without 
an excitation gap. This property persists even if we 
include the bond-alternation, as can be seen in Figs.^ 
and ^ However, it goes without saying that because we 
have such a spin-wave spectrum, we cannot conclude a 
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FIG. 4. The same figure as Fig.^, but with 71 = 72 = —0.9. 

gapless spectrum of the system via the lowest-order spin- 
wave analysis. This is due to the fact that higher order 
quantum corrections play a crucial role whether or not 
the system is indeed gapless. What we would like to 
emphasize in this section is that the lowest mode has a 
linear spectrum, for which we are ensured to use non- 
linear sigma model techniques in order to investigate the 
characteristic properties of the models. 

IV. NON-LINEAR SIGMA MODEL APPROACH 

As was shown in the previous section, the spin excita- 
tions in both models A and B are characterized by the 
linear spectrum in the low energy regime. This suggests 
that we can map these models to 0(3) non- linear sigma 
model (NLSM) in order to incorporate the quantum ef- 
fects in a leading order approximation. Although the 
NLSM has been successfully applied to the usual uni- 
form spin chain |^-^, there appear in the present mod- 
els some other massive modes in addition to the lowest 
massless spin excitations, making the problem somewhat 
complicated. Therefore, it is not straightforward how to 
incorporate the effects of massive modes into the massless 
mode to obtain the effective NLSM. In this section, we 
propose a way to treat such effects, following the method 



recently developed in ref. for spin-ladder systems. 
The starting action is 

S = Sb + Sh, (4.1) 

where 

1 r 

In the following, we would like to treat the two models 
separately. 

A. Model A 

For the Model A, the spin- wave spectrum has a three- 
band structure when si 7^ S2. Therefore, we introduce 
the sigma-model fields as, 

n(3j + a) = m(3j + a) + {-f'+'^+^aolaii] + a), (4.3) 

where oq is a lattice constant, j = 0, 1, • • • , iV — 1 and 
a = 1,2,3. Note that whether the present models are 
massive or massless is essentially determined by the coef- 
ficient of the topological term, and may not be sensitive 
to how many fields we consider in the effective theory. 
However, it is expected that the introduction of three 
kind of fields /i, I2 and ^3 improves our low-energy effec- 
tive theory, by incorporating the effects due to the mas- 
sive modes to some extent: e.g. a better approximation 
may be obtained for the spin-wave velocity, as is seen 
below. In what follows we assume that the several fields 
thus introduced are sufhciently smooth functions so that 
we can take the continuum limit of them. 

First, let us derive the expression for the Berry phase 
term in the continuum limit, 

3 Af-l 
a=l j=0 

X {uj[n{6j -I- a -f 3)] - uj[n{6j + a)]} 

3 N-l 
a=l ]=0 

X JdTSn{6j + a) ■ [n{6j + a) x drn{&i + a)], 

(4.4) 

where Sa = si,S2,S2 for a — 1,2,3, respectively, and 
5n{&j + a) = n{6j -I- a -I- 3) — n(6j + a). We re- 
place the difference by the differentiation, Sn{6j + a) ^ 
3aodxm{6j + a) + (— )"2aoia(6j + a) + 0{al), where oq is 



4 



the lattice constant. Substituting this formula, we then 
have 

= J d^xm ■ {dxinx drfn) 

a 

{rnxdrm), (4.5) 

where s — si and 

1 / " \ 

(4.6) 




with a defined in eq.(3.6) 



Let us next take the continuum limit of interaction 
terms, 



3 



+ 1 



a=l j=0 



X J dr [n{3j + a + 1) - n{3j + a)] . (4.7) 

Here, note that n{3j + a + 1) — n(3j + a) ^ aQdxm{3j + 
a) + ao(-)3j'+'^ [/a+i(3j + a) + ^,(3^ + a)] + 0{al). Sub- 
stituting this expression and neglecting oscillating terms, 
we have 



Sh 



JaoS2 



h], (4.8) 



where K = (2a + l)/3, = (0,0,0) and 

2a a a 

L = - \ a 1 + a 1 

a 1 1 + a 



(4.9) 



Here we have introduced g for later convenience (Model 
B). 

Integrating the field I, we end up with the effective 
Lagrangian density for m. 



where 



V 



OTT 



1 1 



(4.10) 



g 



V — JaoS2\ 



iK-lgtL-^g 



(4.11) 



provided that L is a symmetric matrix. For the Model 
A, by the use of the following inverse matrix of L defined 



in eq.(4.9) 




a + 2 —a —a 
—a a + 2 a — 2 
-a a-2 a+2 



(4.12) 



we have the following formula, 

9 = 27risi, 

2 3 



S2 ^(2a + l)(a2-2a + 4)' 
V — 2JaoS2i 



2a + 1 
■~2a + 



(4.13) 



It should be noted that in the case of si = S2 = s, the 
above formulae reduce to those of the uniform spin-s spin 
chain, 9 — 21:13, g = 2/s and v = 2JaQS. 



It is interesting to compare the formula ( 4.13 ) as those 
in Appendix A, which was previously obtained by assum- 
ing a single Z-field (i.e., neglecting the effect of massive 
modes) The coefficient of the topological term is 

the same, 9 = 9, which makes the above formula reliable. 
The difference of the bulk quantities appears in the fac- 
tors depending on a. From Table ^, we found v closer to 



Sl 


1/2 


1 


3/2 


2 


5/2 






1.000 


1.061 


0.878 


0.707 


0.582 


S2 = 1/2 


V 


1 


1.113 


1 


0.866 


0.761 




V 


1 


1.25 


1.4 


1.5 


1.571 




^sw 


1.488 


2.000 


2.179 


2.121 


1.952 


S2 = 1 


V 


1.569 


2 


2.219 


2.236 


2.138 




V 


1.6 


2 


2.286 


2.5 


2.667 



TABLE I. Spin- wave velocity for the Model A. The veloc- 
ity denoted by fsw is calculated by the spin-wave dispersion 
rel ation u; ~ iigwP for small p. It is compared with v given by 
eg. ( 4. 131 ) in sec. IV. For reference, we also show v given by 
eq.(Al[) for which the effects of massive modes are neglected. 

the value Vsw We can thus see that the approximation 
is improved by the introduction of three /-fields rather 
than a single Z-field. 

Based on the above observations, we now arrive at the 
following conclusion for the Model A. The topological 
term is controlled only by the spin si, and this char- 
acteristic property is independent of the approximations 
we have used. Namely, if si = integer (half-integer), the 
model is to be massive (massless). This is indeed con- 
sistent with the result by Lieb-Shultz-Mattis theorem in 
sec. II. 



5 



B. Model B 

In the spin- wave analysis, we have found a two-band 
structure in the spectrum for the Model B, and so we 
here introduce the fields, 



n{j)^m{j) + {-y+'^aoX 



for 



J = 1,2 
J =3,4 



(mod 4) . 
(4.14) 



Manipulations for passing to the continuum limit are per- 
formed in parallel to the Model A. First, note that the 
continuum limit of the Berry phase terms takes the same 



form as eq.(4.5), but with 

Sl + S2 



f 



(4.15) 



Here, subscript a in eq.(4^) runs from 1 to 2, because 
we have two modes now. Next, inte ract ion terms are 
also calculated in the same form as eq.(4^) for which the 
corresponding parameters are given as 



2Ta 



g = 



Ta — T2 
, 1 / 2ria2 + 
^ - 2 I Ta 



Ta 



Ta 
2T2 + Ta 



(4.16) 



where = 1 - 7, and T = 1 + (71 + 72)72. 

It is straightforward to derive the explicit formulae for 
general 71 and 72, by eq.( 4.11 ). However, their appear- 
ance is rather complicated generally, so that we here write 
down them for the simple case 71 = 72 = 7. We thus have 



= ni{si + S2) 



4(1 + 7)siS2 



1 



9 



(si + S2y +7(si - 52)^ ' 
(si + S2)^ +7(si - 82)" 



S1S2 ^(1 + 7)[(si + S2)2 + -f{si - 6siS2 + S2)] ' 

V = 4:JaoSiS2\ 



(l-7)(l-7^) 



(si + .S2)2 + 7(si - 6siS2 + sj) ' 



(4.17) 



Now we wish to discuss the above formulae in several 
limiting cases. First, by setting si — S2 — s, they are 
reduced to 6* = 27ris(l 7), g = 2/s x (1 - 7^)^^/^ 
and V = 2Jaos(l — 7^)^^^, i-e, the usual formula for the 
uniform chain with bond-alternation. Next set 7 = for 
arbitrary si and 82- Then we have 



— Airi- 



S1S2 



Sl + S2 
Sl + S2 



S1S2 



S1S2 

Sl + S2 ' 



for 7 = 0, 



(4.18) 



We predict from this expression that 7 = the sys- 
tem is massive in general. For example, 9 = ini/S in 
the case si = 1/2 and S2 — 1- Note that the bulk quan- 
tities g and v are different from g and v, respectively, 
defined by eq. ( [A^ ) which were previously derived ||l5| by 
neglecting massive modes from the beginning. On the 
other hand the coefficient of the topological term is the 
same in both cases, 9 — 9. By observing that v repro- 



Sl 


1/2 


1 


3/2 


2 


5/2 






1.000 


1.333 


1.500 


1.600 


1.666 


S2 = 1/2 


V 


1 


1.333 


1.5 


1.6 


1.667 




V 


1 


1.491 


1.936 


2.332 


2.687 




^sw 


1.333 


2.000 


2.400 


2.666 


2.857 


S2 = 1 


V 


1.333 


2 


2.4 


2.667 


2.857 




V 


1.491 


2 


2.500 


2.981 


3.440 



TABLE II. Spin-wave velocity for the Mo del B . Parameters 
are same as Table ^ but with v given by eq.(4.17), and v given 
by eq.(0). 



duces Vsw exactly in Table O, we naturally expect that 



the assumption (4.14) in our approach works rather well, 
making the above formula more reliable than the previ- 
ous treatment with a single-/ field approximation. 

Next, let us discuss the effects of the bond-alternation. 
First, we recall that for the usual uniform spin-s chain 
(see the above formula for si = S2 case), 6 = 0^ Anis 
if 7 varies 7 = — 1 — s- 1. Namely, the model realizes 
2s massless points during this process [ p2[ . However, it 
is difficult in general to predict the exact value of the 
bond-alternation parameter 7 which brings about mass- 
less phases. For example, when s = 1, it is predicted 
by Affleck and Haldane |2^] that a massless phase occurs 
at 7 = 1/2, while the numerical calculation suggests a 
different value 7 ~ 0.25 ||2^]. This is because for systems 
with explicitly broken translational symmetry, there are 
no strong restrictions by this symmetry and the coeffi- 
cient of the topological term can take arbitrary values. 
This makes such values less reliable. Nevertheless it is 
believed that the number of massless points predicted by 
sigma model techniques itself should be correct, which 
has been indeed realized for the s = 1 case by numerical 
calculations ||23|| . 

Keeping this in mind, we would like to observe how 
the topological term behaves for for si ^ S2 case. From 
eq.( 4.17| ), we can see that 9 is an increasing function of 
7) 



= ^ 7ri(si + S2) 



4siS2 



(4.19) 



according to 7 = — 1 — > 1. For example, in Fig.|^, we 
explicitly show 9 (see the curve (a) in the figure) as a 
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FIG. 5. (a) 9 and (b) 9 in unit vri as functions of 7 for 
Si — 1/2 and S2 = 1. 

function of 7 in the case si = 1/2 and §2 = 1. We can 
see that at 7 ^ —0.273, 6* becomes 1 and massless phase 
is expected. It is interesting to compare this value with 
the numerical one 7 ^ —0.13 in ref. |jl^. This difference 
is not serious, as has been already discussed. What we 
are interested in is mainly how many times we expect 
massless phases in varying the bond-alternation parame- 
ter. Indeed, the conclusion of Tonegawa et al |1^ is that 
there occurs only one massless phase in the whole anti- 
ferromagnetic region of 7, which is consistent with our 
results. 

As we have now two kind of spins, it is interesting to 
study an extreme case where the difference between the 
two spins is large. Set, e.g, S2 ^ 00, and we have 9 — s- 
Airisi (and g — > 0) near 7 = 1. In this limit, therefore, we 
have 2s 1 massless points during the process 7 = — 1 ^ 1. 
Namely, it turns out that smaller spin solely controls the 
number of the massless points. 

Finally, we would like to compare the present result 
with the previous one (A2) which was derived by ne- 
glecting the effects of the massive spin modes. It turns 
out from ( |A^ ) that 9 = 0^ 27ri(si -I-S2) for 7 = -1 ^ 1. 
Therefore, in the case that si and S2 are the same type 
of spins (integers or half-integers) the model has si + S2 
massless points, while it may have si + S2 — massless 
points in the case where the spins are of different type 
(the 7=1 point cannot be classified as a massless point 
because the system is separated into disconnected dimers 
there). In Fig.^ we present 9 for reference, from which 
we see that behavior near 7 = 1 is different from that of 
9. However, the value of 7 which gives 6* = 1 is similar to 
that for 9 in this case. In general, 6* < 61 at 7 = 1, where 
equality holds only if si = S2. This discrepancy becomes 
more serious if the difference between the two spins si 
and S2 becomes larger. 

Then the question is which formula is more reliable. In 
table [II, we compare the spin- wave velocities calculated 
by means of the above two methods in case of the finite 
bond-alternation, corresponding to Figs.^ and IJ. It is 
clearly seen from this table as well as Table || that not 





1/2 


1 


3/2 


2 


5/2 








0.436 


0.335 


0.279 


0.255 


0.242 


7 


= 0.9 


V 


0.4JD 


O.oot) 


0.2/9 


U.ZOO 


U.242 






V 


0.436 


0.722 


1.089 


1.504 


1.952 






'^sw 


0.436 


0.614 


0.748 


0.859 


0.955 


7 


= -0.9 


V 


0.436 


0.614 


0.748 


0.859 


0.955 






V 


0.436 


0.591 


0.689 


0.768 


0.837 




TABLE III. 


Spin- wave 


velocity 


for the 


Model 


B with 



bond-alternation. The spin S2 is fixed as S2 
parameters are same as Table ^ 



1/2, and other 



only for 7 = but also for rather large 7, w in ( 4.17 ) 
exactly reproduces the spin-wave velocity Vsvi calculated 
in the spin-wave analysis in the previous section. On 
the contrary, the velocity v considerably deviates from 
Vsw Notice especially the difference between Figs. || and 
^. These figures tell us that the hybridization between 
two spin modes in Fig.^ is much larger than that in Figj^. 
Therefore, if we naively neglect the massive mode, it may 
cause more serious errors for the former case. One can 
indeed see this fact in the above table: Namely, w's for 
7 = 0.9 is worse in general than those for 7 = —0.9. 
Based on these observations, we believe that the present 



formulae eq.(4.17), which incorporate the effect of mas- 
sive modes, provide more reliable results. 

To conclude, we may predict how many massless 
phases we encounter when we vary 7 from —1 to 1: 
Namely it is given by the number of odd-integers included 
in the region between and 4siS2(si -f S2)/{s\ + s^). It 
may be an interesting issue to investigate this conjecture 
directly by numerical calculations. The simplest case is. 



CD 2 





(bh/ 











-0.5 



0.5 



FIG. 6. (a) 9 and (b) 9 in unit ixi as functions of 7 for 
Sl = 1/2 and S2 = 3/2. 

for example, si — 1/2 and S2 = 3/2. The present formula 
{9) predicts only one massless phase, while the previous 
formula (9) predicts 2 massless phases, as can be seen 
from Fig.pl 
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V. SUMMARY 

We have investigated two kind of alternating spin 
chains, both of which have singlet ground state. They 
are composed of a periodic array of two kind of spins. To 
analyze low-energy properties of these models, we have 
first studied spin-wave spectrum in detail. Introducing 
quantum fluctuations semiclassically, we have converted 
the low-energy spin mode to the 0(3) non- linear sigma 
model. Analyzing the topological term in the sigma 
model, we have then clarified how the massless or mas- 
sive excitations are developed reflecting the topological 
nature of the alternating spin system. Up to now, the 
alternating spin chains found experimentally have ferri- 
magnetic ground state. We think that the systems simi- 
lar to those proposed here with singlet ground state could 
be fabricated experimentally, they would provide a new 
paradigm of the quantum spin systems, for which one 
may systematically observe the competition of the mas- 
sive and massless phases. 
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APPENDIX A: FORMULAE WITHOUT 
INTRODUCING MASSIVE MODES 

In this appendix, in comparison with the formulae 
given in the text, we write down the expression with- 
out taking into account the effects of massive modes. We 
here quote the results obtained in the previous paper 
, in which a preliminary treatment was given for low- 
energy properties of the model. There, in order to map 
the model to the non-linear sigma model, only a single 
Mield was introduced, which means that the effects due 
to massive modes are neglected. One can find the detail 
of the derivation is given in the paper . The formulae 
obtained are 

6 = 27risi, 
9 



S2a + 2' 

2a + 1 



V — 2JaQS2 



a + 1 



for Model A, and 



27riscff(l -I- 7off), 



cffVl - 7, 



2 ' 
off 



V = 2JaoSoff-y/l - 7ofT: 
for Model B with 71 = 72 = 7, where 

Sl + S2 



(A2) 



Scff = 



7off 



(51 + 52)^7 " jsi - 82)^ 



(A3) 
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